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Abstract 

We consider the Langevin dynamics of a Gaussian test polymer chain 
coupled with a surrounding matrix which can undergo the glass tran- 
sition. The Martin-Siggia-Rose generating functional method and the 
nonpertubative Hartree approximation are used to derive the general- 
ized Rouse equation for the test chain. It is shown that the interaction 
of the test chain with the surrounding matrix renormalizes the bare 
friction and the spring constants of the test chain in such a way that 
the memory function as well as the bending dependent elastic modulus 
appear. We find that below the glass transition temperature Tq of the 
matrix the Rouse modes of the test chain can be frozen and moreover 
the freezing temperatures (or the ergodicity-nonergodicity transition 
temperature) T c (p) depends from the Rouse mode index p. 
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I. INTRODUCTION 



It is wellknown that for relatively short polymer chains the standard Rouse 
model can describe the dynamics of a melt reasonably well On the contrary, 

for chain length N exceeding a critical length, the entanglement length N e , the 
behavior is usually described by the reptation model Here we restrict ourselves 
to chain lengths iV < N e , i.e. the entangled polymer dynamics will be beyond of 
our consideration. 

The reason why in a dense melt the Rouse model provides so well dynamical 
description for short chains is connected with a screening of the long-range hydro- 
dynamic as well as the excluded volume interactions. As a result the fluctuations 
of the chain variables are Gaussian. But there are further essential questions: How 
does the bare monomeric friction coefficient £o an d the entropic elastic modulus e 
(which are simple input parameters of the standard Rouse model) change due to the 
interactions of the test chain and the surrounding matix? Why does such a simple 
model work so well for describing short chain melts? Obviously, the corresponding 
answers cannot be given by the Rouse model, which describes only the dynamics of 
connected Gaussian springs without further interactions. 

On the other hand, at relatively low temperatures close to the glass transition 
of the surrounding matrix the deviations from the standard Rouse behavior will 
be definitely more pronounced. For example, Monte Carlo (MC) studies of the 
bond fluctuation model at low temperatures (but still above the temperature region 
where possibly the glass transition mode coupling theory [[| applies) show that the 
Rouse modes remain well-defined eigenmodes of the polymer chains and the chains 
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retain their Gaussian properties . Nevertheless, the relaxation of the Rouse modes 
displays a stretched exponential behavior rather than a pure exponential. It could 
even be expected that at temperatures below the glass transition temperature of 
the matrix Tq the Rouse modes are frozen out. In these temperature regimes the 
interactions between monomers take a significant role and determine the physical 
picture of the dynamics as will be shown below. 

The generalized Rouse equation (GRE), which can be used for the investigation 
of the problems mentioned above, has been derived by using projection formalism 
methods and mode coupling approximations (MCA) As a result of projection 

operator formalism the time evolution of the test chain is expressed in terms of a 
frequency matrix, which is local in time, and a memory function contribution due 
to the inter-chain forces exerted on the test chain segments. With the assumption 
that the frequency matrix term has the same form as in the standard Rouse model 
(linear elasticity with the entropic modulus e = 3kf,T/l 2 ) all influence of the matrix 
chains reduce to the memory function contribution 0-0]. 

The projection operator methods appears to be exact but rather formal, and to 
derive explicit results further approximations have to be made, which can be hardly 
controlled often. Therefore it is instructive to use another alternative theoretical 
method to derive the GRE. Recently, a non-pertubative variational method which is 
equivalent to a selfconsistent Hartree approximation was used for the investigation 
of the dynamics of manifolds || and sine-Gordon model in a random media. 
As a starting point the authors employed the standard Martin-Siggia-Rose (MSR) 
functional integral technique |IO,lT|. Here we follow this approach to derive a GRE 



and study the dynamics of a test polymer chain in a glass forming matrix. 



The paper is organized as follows. In section 2, we give a general MSR- functional 
integral formulation for a test chain in a polymer (or non-polymer) matrix. Under 
the assumption that the fluctuations of the test chain are Gaussian the Hartree- 
type approximation is applied and a GRE is finally derived. The case when the 
fluctuation dissipation theorem (FDT) and the time homogenity are violated is also 
shortly considered. In section 3 on the basis of the GRE some static and dynamical 
properties of the test chain are discussed. In particular the theory of the test chain 
ergodicity breaking (freezing) in a glassy matrix is formulated. Section 4 gives some 
summary and general discussion. The appendices are devoted to some technical 
details of the Hartree-type approximation. 

II. GENERALIZED ROUSE EQUATION (GRE) 

A. MSR-functional integral approach 

Let us consider a polymer test chain with configurations characterized by the 
vector function R(s, t) with s numerating the segments of the chain, < s < N, 
and time t. The test polymer chain moves in the melt of the other polymers (matrix) 
which positions in space are specified by the vector functions r^(s, i), where the 
index p = 1,2, ...,M numerates the different chains of the matrix. The test chain 
is expected to have Gaussian statistics due to the screening of the self-interactions 
in a melt 0. We consider the simultaneous dynamical evolution of the R(s,t) and 
r( p )(s,t) variables assuming that the interaction between matrix and test chain is 
weak. 

The Langevin equations for the full set of variables {R(s, t), rW(s, t), . . . , r^ M ^(s, t)} 
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has the form 



Co^M) - s-R 3 (s,t) + — -- £ / ds'V (R(s, t) - r^s't) 



d 



M rN 



5Rj(s,t) f^Jo 



fi(s,t) 







2 



(p), 



5 ^ 



ds 2 



pl\ 

£ / ^Y(rW(s,t)-rH( s ^) 

™ 1 ** 



+ 



5 ^ 



AVJ - ply 

E / dsV(r^-R(^ = 

71=1 J0 



(2) 



where £o denotes the bare friction coefficient, e = 3T/1 2 the bare elastic modulus with 
the length of a Kuhn segment denoted by I, V(- ■ ■) and V(- ■ •) are the interaction 
energies of test chain-matrix and matrix-matrix respectively, and fj(s, t), fj(s, t) are 
the random forces with the correlator 

a(M)£(s',0)> = (/i(M)/i(*',0) = 2r£oM(*-*W-0- (3) 



After using the standard MSR-functional integral representation |T0[ for the system 
([l]-^), the generating functional (GF) takes the form 
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where the dots represents some source fields which will be specified later and Ein- 
stein's summation convention for repeated indices is used. In GF (^) the MSR-action 
of the free test chain is given by 
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As we will realize later the explicit form of the full action of the medium 
A\ r^(s, t), r ^'{s, t) plays no role. In principle it could have any form and in 
particular, for a polymer matrix, the following one 

d i n \ , d 2 
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p =i Jo J 

In order to obtain an equation of motion for the test chain one should integrate over 
the matrix variables r^(s, t) first. For this end it is reasonable to represent GF (Q) 



as 



£{•••} = / DR^s^DR^s.t) 
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where the influence functional H is given by 
s R, R 



-In / H Dr^\s,t)DT^\s,t) 
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In the spirit of the mode coupling approximation (MCA) the force between 
the test chain and the matrix should be expressed as a bilinear product of the two 



subsystems densities. In order to assure this we expand the influence functional 



with respect to the forces Fj = — V jV between the test chain and the matrix up to 
the second order. This leads to 
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where the matrix density 

i Jo 



M . N 



(10) 



and the response field density 

M „ N 

n,(r,t) = E/ dstrf\s,t)5(v-T^(s,t)) (11) 
p=i ^ 

were introduced and (■ ■ -) x denotes cumulant averaging over the full MSR-action 
Ai [r, f] of the matrix. In eq. @ the term (f <^ t) is the same like the previous one 
but with permutated time arguments. The terms which are linear with respect to Fj 
vanishes because of the homogenity of the system. In the Appendix A we show that 



because of causality the correlator (H(r, t)ILj(r', t')), equals zero ||10|-|T2|1. Taking 
this into account and performing the spatial Fourier transformation the expression 
for GF (^) takes the form 
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£{• •-} = J DRj(s, t)DRj(s, t) exp j - A [R(s, t), R(s, t) 

+ \J dsds ' J dtcMiRjfat) J -^^kjhlVik)] 2 S(k,t - 1') 

x exp {ik [R(s, i) - R(s', t')]}iRi{s', f) 
+ J ds ds J dtdt' iRj(s,t) J -^^kjki\V(k)\ 2 Pi(k,t - t') 

xexp{ik[R(s,t)-R(s',t / )]}| ( 12 ) 
where the correlation function 

5(k,t) = (p(k,t) /) (-k,0)) 1 (13) 

and the response function 

p,(k J t) = <p(k,t)n,(-k,o)) 1 (14) 

of the matrix are naturally defined. Going beyond the LRT-approximation would 
bring us multi-point correlation and response functions. 

We should stress that in contrast to the matrix with a quenched disorder which was 
considered in |||§ in our case the matrix has its own intrinsic dynamical evolution 
which is considered as given. For example, for the glass forming matrix, which is 
our prime interest here, the correlation and response functions are assumed to be 
governed by the Gotze mode-coupling equations ||. 

B. The Hartree approximation 

The Hartree approximation (which is actually equivalent to the Feynman varia- 
tional principle) was recently used for the replica field theory of random manifolds 



[ |T3|j as well as for the dynamics of manifolds || and sine-Gordon model || in a 
random media. 
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In the Hartree approximation the real MSR-action is replaced by a Gaussian 
action in such a way that all terms which include more than two fields Rj(s, t) or/and 
Rj(s, t) are written in all possible ways as products of pairs of Rj(s, t) or/and Rj(s, t), 
coupled to selfconsistent averages of the remaining fields. As a result the Hartree- 
action is a Gaussian functional with coefficients, which could be represented in 
terms of correlation and response functions. After these straightforward calculations 
(details can be found in the Appendix B) the GF (|I2D takes the form 

Z{---} = J DR j (s,t)DR j (s,t)expS^ - A [R(s, t), R(s, t) 

pN poo c-t 

+ / dsds' dt dt' iR j (s,t)R j (s',t')X(s,s';t,t') 

J0 J—oo J—oo 



pN poo a pt 

— / dsds' dtiRj(s,t)Rj(s,t) / dt f \(s, s';t,t') 

JO J—oo J—oo 

| pN p°° ft ^ a "1 

+ -/ dsds' dt d£iRj{s,t)iRj{^,l!)x{s,s'',t,i!) } (15) 

2 JO J-oo J-oo 



where 
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In eq. ([16| , |T7D the response function 

G{s,s';t,t') = (i&{s',t')'R(s,t) 

the density correlator 



F(k; s, s'; t, t') = exp j ~ Q(s, s'; t, t') \ (19) 



with 



Q(s,s';t,t') = (R(s,t)R{s,t)) - (R(s, t)R(s', t')) 
= C(s,s;t,t) - C(s,s';t,t') 
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and the longitudinal part of the matrix response function 



P(k;t,t') = ikjPjfct,?) 



(21) 



are defined. The pointed brackets denote the selfconsistent averaging with the 
Hartree-type GF flOf). 

Up to now we considered the general off-equilibrium dynamics with the only restric- 
tion of causality [flO|— p72|| . We now assume that for very large time moments t and 
t', where the difference t — f is finite so that — > 0, time homogenity and the 
fluctuation-dissipation theorem (FDT) holds. This implies 
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where (3 = 1/T. By using this in eq. flT5| ) and after integration by parts in the 
integrals over t' the GF in Hartree approximation takes the form 
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where the subscript "st" indicates the static correlation functions. This generating 
functional immediately leads to the following generalized Rouse equation (GRE) 

d r N r l d 

£o^i*i(M) + I ds' J ^{s.s'-t-t^—R^s'^) 

r N 
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where the memory function 
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and the effective elastic susceptibility 
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are defined. The correlation function of the random force T« is given by 



(26) 



(27) 



(Fi(a, t)Fj(s', t')) = 2T5ij k 6(s - s')5(t - t') + 6(t - t')T(s, s'; t - t') 



:28i 



As a result we have obtained basically the same GRE as in the papers but 
with one additional elastic term. This term (see the 2-nd term in eq. (|27|)) is mainly 
inversely proportional to the temperature and is, in contrast to the first term, of 
an energetic nature. The two factors of kV(k) quantify the forces exerted by a 
pair of surrounding segments on the test chain segments s and s', whereas the 
S st (k) and F st (k; s, s') factors quantify the static correlations between the segments 
of surrounding and test chain segments, respectively. In |5|-[7| only the entropic 
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elastic part was taken into account. The memory function (|26|) has the same form 
as in [§-0] and the relationship fl28|) is assured as soon as the FDT fl22|) and (p3|) is 
fullfilled. 



C. Generalized Rouse equations for the off-equilibrium dynamics 



In this subsection we give GRE's for the more general case when the time ho- 



mogenity (stationarity) and the FDT do not hold [13]. 

By employing the standard way || one can derive two coupled equations of motion 
for correlators C(s, s'; £, t') and response functions G(s, s'; t, if) 
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with the initial conditions 



£ G(s,s';t = t' + + ) =<J(s-s') 
G{s,s';t = t')=0, 



t < t' 



and 



C{s,s'\t = if) = (R{s,t)R{s',t)) 



(30) 



(31) 



(32) 



In the stationary case all correlators and response functions in eq. ( |2l)| - |3"2"|) only 
depend from the differences of time moments, t — t'. If we assume again that FDT 
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2D and (|23]) holds, then from eq. ( p0[) after performing the integrations by parts 
in the integrals over t") one arrive at the GRE for t > 

Z -C(s,s';t) + j Jt> j o ds"T(s,s";t-t')—C(s",s';t') 

- / ds"tt(s,s")C(s",s';t) = 0. (33) 
Jo 



Of course, eq. (P3f) could be obtained immediately from eq. (EH) by multiplying both 



sides with Rj(s', 0), averaging and taking into account that because of causality 



(J r (s,t)R(s',0)) = at t > 0. We will use the GRE eq.©, where the functions T 



and f2 are given by eqs. ( Pq , |27| ), in the next section for the investigation of the test 



chain ergodicity breaking (freezing). 

III. SOME STATICAL AND DYNAMICAL PROPERTIES OF THE TEST 

CHAIN 

The new features of the GRE (^) relative to the standard Rouse equation are 
that it contains the integral convolution with respect to the s— variable in the fric- 
tional term as well as in the elastic term. The frictional term is also non-local in 
time. All these things together should change the statical and dynamical behaviour 
of the Gaussian test chain in comparison to the ideal chain. 

We also should stress that the GRE is substantially nonlinear because the memory 
function (|26| ) depends from the test chain correlator C(s,s';t) in such a way that 
a positive feedback obviously exists. That is the reason why one could expect that 



eq. (|33|) shows an ergodicity breaking in the spirit of Gotze's glass transition theory 
I- 

As usual it is convenient to introduce the standard Rouse mode variables 
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X(p, t) = ± j" dsR(s, t) cos (^f) (34) 



with the inverse transformation 

R(5, t) = X(p = 0, t) + 2 ]T X(p, t) cos . (35) 

In general one also needs a 2- dimensional Rouse transformation 

r (P>9;*) = ^ / rfs'T(s',s";t)cos ^^-J cos ^^-J (36) 

where functions like T(s',s") should be treated like N x iV-matrices. For example 
the density correlator (|1£]) should be considered as an exponential function from a 
iV x iV-matrix Q(s , s ) and the series expansion holds : 
P 1 /P\ 2 fiV 

F( s , s ') = i-yQ( s , s ) + -fyj y o d s g( S)S )Q(s,s) 

1 /h 2 \ 3 rN r N 

-31(3-) 7 rfs 7 rfs )Q(a + (37) 

We also assume that matrices in the Rouse mode representation are nearby diagonal 

r(p,g) = <y p , ff r(p) (38) 

Qip,q) = S m Qip) (39) 

n( P ,q) = 6 M Q( P ) (40) 

for any p and q not equal zero |J . 

Then as a result of Rouse mode transformation the GRE for the Rouse mode time 
correlation function, C(p,t) = (X(p, t)X(p, 0)), takes the form (for p ^ 0) 

Zofc( P: t) + J*dt'T( P ,t-t')-^C(p,t') + n(p)C(p,t) = (41) 

where 
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T(p 7 t) = (3 J -^k 2 \V(k)\ 2 lexp 



(2vr 



-NC(p,t) 



iU(k,t) 



and 



^iVWI ^k) F st (k;p)-F st (k;p = 0) (43) 



For p = the GRE describes the dynamics of the centre of mass 



1 r N 

R c . m (t )=X(p = 0,t ) = — / dsRs.t) 

A/ Jo 



and has the following form 



£ Q |R, m (t) + /? jT <tf / ^ fc 2 ] | 2 F(k; p = 0, 9 = 0, t - t') 



d 3 k 
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with 



1 /" w 

(f cm) jit) = J^J dsTj{s,t) 



and 



1 fJV /• 

F(k;, = 0, g = 0;t) = -/ o <W 



(42) 



(44) 



(45) 



(46) 



(47) 



As a result all Rouse mode variables relax independently. The conclusion that the 
Rouse modes are still "good eigenmodes" even in the melt is supported by Monte- 
Carlo @] and molecular-dynamic Jl5j simulations. 



For cases where the assumption of diagonality [38,39,^0| cannot be justified, the 
Rouse modes do not decouple and one have to go back to eq. (R3). In the Rouse 



mode representation it reads as 



^jC( P ,t) + jfdf JdqT(p,q;t-t')^C(q,t')+ J dq n(p, q)C(q,t) = (48) 
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A. Static properties 

As we have already discussed in sec. II. B the interaction with the surrounding 
segments renormalizes the elastic properties of the Rouse chain so that the test 
chains elastic susceptibility is given by eq. ([431). The additional elastic term in GRE 
leads to the renormalized static normal modes correlator 

CM = Jim (49) 

Explicit evaluation of the Q(p) can be done if we use for the static correlator F st (k; p) 
the standard Rouse expression 

> Vs \os(^). (50) 



1 r N ( 

F st (k;p) = — J dsexp I - 



6 J V N J ' 
Then the calculation yields for the two limiting cases 



n(p) 



where 



(77) [e + £,Cdkk*\V(k)\*S st (k)g(eRl)]+0(p*) : f «1 (50a) 
(f)^+4&/r^fc 4 l^)| 2 ^(A;)(^)(l-e-^) : £- 1 (506) 



2 



2 - (x 2 + 2x + 2)e _x ] , R g = ^r- (52) 



and we have chosen I 1 as a cutting parameter. It is evident from the previous eqs. 
(50a, 506) that at small p 

• the elastic modulus gains an energetic component which, in contrast to the 
entropic part s, increase with the cooling of the system, 

• initially absolutely flexible chains acquires a stiffness because of terms of order 
p 4 and higher. 
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At large p the elastic behaviour reduces to the standard Rouse one, as it is expected. 
In Fig.l is shown the result of a numerical calculation of the static correlator ((49|) . 
The Fourier component of the potential is taken, as it is customary e.g. in the 



theory of neutron scattering fll6| , in the form of a pseudo potential approximation, 
V(k) = 7<t 3 , where 7 and a have dimensions of molecular energy and distance, 
respectively. The static structure factor S' st (k) is chosen in the form of the Percus- 



Yevick's simple liquid model |T7[]. One can see that for iV = 500 the small Rouse 



mode index limit (50a) starts at < 3 ■ 10~ 3 whereas the opposit limit (506) is 
fullfilled at > 10 _1 . Because the correlator C st {p) depends mainly from p/N, 
for relatively short test chains the high mode index limit (506) is shifted into the 
window of calculations (see Fig.l for N=20). 

At least qualitatively this deviation from the standard Rouse behaviour have been 
seen by Kremer and Grest in their MD-simulations (see Fig. 3 in ||15||). 



B. The test chain ergodicity breaking in a glass forming 

matrix 

First we consider the case p 7^ 0. In the nonergodic state the Rouse mode 
correlation functions can be represented as 

*(p,f) = = *«*(P.*) +9(P) (53) 

Cst(p) 

where the non-ergodicity parameter 

g(p)=)imV(p,t) (54) 

t— >oo 

was introduced and ^ reg (p,t — > 00) = 0. 
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For the correlation function of the glassy matrix we can use the standard result 
of the glass transition theory || 

^•^fM^ r(k) + ' ,(k)A ' /2 (T)° < 55 > 

where the proximity parameter A = (Tq — T)/Tq is defined and Tq is the temper- 
ature of the matrix ergodicity breaking (Gotze temperature). In eq. ( p5| ) / c (k) is 
the non-ergodicity parameter of the matrix, oc A _1 / 2a is the characteristic time 
scale, a is the characteristic exponent, < a < 1/2 and h(k) is some amplitude. 
In order to derive the equation for g(p) let us take the limit t — > oo in eq. fl4l[) 
keeping in mind the definitions (|54l) and (|55|). Very close to the test chain ergodicity 
breaking temperature T c (p), g{j>) goes to zero (A-type transition 0) and we can 
expand the exponential function in eq. fl42|) up to the first order with respect to 
g(p). The solution of the resulting equation has the simple form 

() = 1 6^(p) 2 

9{p> j^k*\v(kws st (k)fc(k) 1 } 

The critical temperature T c (p) is determined by the equation 

g(p,T = T c ) = (57) 



The numerical solution of eq. (p7j ) is given in Fig.2. It is obviously that if the 
entropic part of fl(p) dominates, the critical temperature is given by 

T c (p) oc (—) . (58) 

Fig.2 really shows that this law ([58]) is well satisfied due to the fact that the critical 
temperatures T c (p) are quite high. But for low temperatures the energetic con- 
tribution in Q(p) is enhanced which leads to a deviation from this simple (N/p) 2 - 
dependence. 
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Now we consider the case for p = 0. The equation fl45 ) for the velocity of the center 
of mass 

V(t) EE ^Rc.m(0 (59) 

leads to the equation for the velocity correlator 

£o (vjitMO)) + f dt'Yit - (^(0^(0)) = ((/c m ),(*k(0)> (60) 



o 



where 



r(t) = /?/ ^A; 2 |\/(A;)| 2 F(k;p = g = 0;t)5(k,t) (61) 

Because of causality the correlator on the r.h.s. of eq. (|60[) has the form 

, t > 

((/cm)i(*)«i(0)> = <( (62) 
^ , t = 



where, as it comes from eq. ([45] 



Vi(0) = y {f c . m )i (0) (63) 



Taking into account the definition of (f c .m)i (t) and eq. (p8|) this yields to the corre- 
lator 

1 r N r N ,, 
(ifcmUtMO)) =2T5 lJ 5(t)— I t^jf ds^s'-s) 

2T 

= -^M(f) ( 64 ) 

Because of the causality property (61) only the ^-functional term on the r.h.s. of 
eq. fl28|) contributes to the correlator (@). Therefore the resulting equation for the 
self-diffusion coefficient 

0=ijHdt<v(t)v(O)> (65) 
19 



takes the form 



D 



T 



N[to + JS°der(t)] 



(66) 



which was obtained before in |5|.[|. 

One can calculate the second term in the denominator of eq. ( |66"D selfconsistently. 
Because now the relevant times t ^> r rouse the approximation 



Q(s', s ;t) = 6Dt + I \s' — s | + const. 



could be used in eq. (flop. Then the density correlator ([47]) is given by 



(67) 



F(k; V = q = 0; t) = m 1 exp (~k 2 Dt) (6* 

12 ™ 

With the use of eqs. (j68|) , (|55|) and eq. fl6ll) in the limit .D — > eq. ( |S6|) becomes 



T 



iV 



til f d3fc |V(fc)PS Jt (k) f /.N 



(69) 



where the denominator is given by static properties only. Similar statements have 



been suggested already in [18 , [19| The solution of eq. (|S9|) has the simple form 



N 

r> = ( i - T2 x 



(70) 



where 



X 



d 3 k \V(k)\ 2 S st (k) 



(2tt)3 ^i + iV 



(71) 



Finally, the temperature of the ergodicity breaking (localization) for the mode p = 
of the test chain is 



T c (p = 0) = (N X ) 1/2 



(72) 
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Fig. 3 shows the results of numerical calculations of T c (p = 0) and T c (p = 1) as 
functions of N. One can see that in the reasonable range of parameters T c (p = 0) > 
T c (p =1). As a result one can say that on cooling of a test chain in a glassy matrix 
the mode p = is the first to be freezed. On the subsequent cooling the modes 
p = 1,2, ... ,N are freezed successively, 

T G > T c (p = 0) > T c (p = 1) > T c (p = 2)>...T c (p = N). (73) 

It is apparent that the system studied here is a nontrivial polymeric generalization 
of the model introduced by Sjogren [pOfl . This model was used for the investigation 



of the /5-peak in the spectrum of glass forming systems 21 



IV. DISCUSSION 

In this paper we have derived a GRE for a test polymer chain in a polymer (or 
non-polymer) matrix which has its own intrinsic dynamics, e.g. the glassy dynam- 
ics ||. We have used here the MSR-functional integral technique which could be 
considered as an alternative to the projection operator formalism |J. One of the 
difficulties in this formalism is the necessity of dealing with the projected dynamic, 
which is difficult to handle with explicitly. On the contrary in MSR-technique the 
dynamic of slow variables is well defined and several approximations which one have 
to employ could be justified. 

In the interaction of the test chain with the surrounding matrix only two-point 
correlation and response functions are involved. In terms of MCA [[J this obvi- 
ously corresponds to the projection of the generalized forces only onto the bilinear 
variables: product of test chain density and matrix density. 

21 



To handle with the action in the GF of the test chain we used the Hartree-type 
approximation (i.e., equivalent to the Feynman variational principle) |8|,|9|,p!3|| , which 
is reasonable when the fluctuations of the test chain are Gaussian. In the case of a 
polymer melt (high densitiy) this is indeed the case due to the screening effects for 
the excluded volume 

The use of the Hartree-type makes the problem that we deal with analytically 
amenable and results in the GRE's for the case when the FDT holds as well as for 
the case when FDT does not hold. In this paper we have restricted ourselves to the 
first case and have shown that the interaction with the matrix renormalizes not only 
the friction coefficient (which makes the chain non-Markovian) but also the elastic 
modulus (which changes the static correlator). The form of the static correlator for 



the Rouse mode variables is qualitatively supported by MD-simulations flgl . 

As regards the dynamical behaviour, we have shown that the test chain in a glassy 
matrix (with the matrix glass transition temperature Tq) undergoes the ergodicity 
breaking transition at a temperature T c (p) < Tq- The critical temperature T c (p) 
could be parametrized with the Rouse mode index p and is a decreasing function of 
P- 

We have considered only the A-type transition which is assured by the bilinear 
term in the expansion of eq. fl42|) . It seems reasonable that keeping the whole 
exponential function in eq. ([42]) might lead to a B-type transition also. The results 
also essentially would change if the off-diagonal elements in the matrix (^) can not 
be neglected (see eq. (f4~5D). In this case only one ideal transition temperature T c 



would be possible. The general theory of a A-type transition was discussed in 23 



This picture of freezing here should not be mixed with a different one, the un- 
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derlying glass transition by itself (e.g. the glass transition of the matrix at T — Tq). 
According to the present view of this phenomenon ||, the spontaneous arrest of the 
density fluctuations is driven by those of the microscopic lengthscale ko, where ko 
is the wave vector which corresponds to the structure factor's main maximum. The 
freezing of these fluctuations then arrests the others through the mode coupling. 
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APPENDIX A: RESPONSE FIELD DENSITY CORRELATOR 

It is more convenient to handle with the spacial Fourier transformation of this 
correlator 

(n,(k,t)ii,-(-k,o)> 1 (Ai) 

M 

= £ / dsds'(tfi p \s,t)tf i j m \s\0)exp{tk[r^\s,t)-^ m \s',0)]}) i 

(p,m=l) 
oo i M , 

E^E J dada' (^f^ p) ( S ,^)^ff ) ( S ^0)[^krW( S ,^)] a Hkr( m )( S ^0)] 6 ) l 



a,b=0 p,m=l • 



Such multi-point cumulant response functions (MRF) were considered in [12]. The 
causality condition for these functions asserts that the time argument of at least one 
r-variable should be the latest one, otherwise this MRF equals zero. Because of the 
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same reason self-loops of response functions vanish []T0 JTT | - MRF's which consists 



only of f-variables also vanish. 

In the case ( |A1|) all time arguments of the r- variables are equal to the corresponding 
time arguments of f-variables and as a result the MRF in eq. ( pYT| ) vanishes. 



APPENDIX B: DERIVATION OF THE HARTREE-TYPE GF 



In order to calculate the bilinear Hartree action, we follow the way mentioned in 



sec. II. B. With these strategy in mind the 2-nd term in the exponent (12) is evaluated 



as 



1 r N 



2 Jo 

l r N 



dsds' / dtdif iR j (s,t)R j (s',t')I 1 (s, s';t,t')S(k;t,t') 



OG 
OG 



+ -/ dsds' I dtdt' iR j (s',t')R j (s,t)I 2 (s,s';t,t')S{k;t,t') 
2 Jo 



1 r N 



OG 
OO 



+ -/ dsds' I dtdt' iR j (s,t)R j (s,t)I 3 (s,s';t,t')S(k;t,t') 
2 i/o 



l r N 



DO 
OO 



+ -/ dsds' I dtdt' iR j (s',lf)R j {s',t')I 4 (s,s';t,t!)S{k;t,t!) 
2 Jo 



l r N 



OO 
OO 



+ -/ dsds' dtdt' iR j (s,t)R j {s l ,t')I 5 (s,s';t,t')S(k;t,t') 
2 Jo 



1 rN POO 

+ -/ dsds' / dtdt' iR j {s,t)R j (s',t')I 6 (s,s';t,t')S{k;t,t') (Bl) 

2 JO J-oo 



where 



h(s,s';t,t') = 



I 2 (s,s']t,t') 



I 3 (s,s';t,t') 
h(s,s';t,t') 



S 



d 3 k 



SR^s'^J (2tt) 3 
5 r d 3 k 



kjklVik^expiik [R(s,t) - R(s',t'))}iRi(s',t' 



5Rj(s,t)J (2tt) 3 



kjki\V(k)\ 2 exp {ik [R(s, t) - R(s, t' )]} iRi(s, t) 



5 



d 3 k 



x 5Rj(8,t)J (2tt) 3 
d 3 k 



5 



5Rj{s',t')J (2tt) 3 
d 3 k 



kjklVik)] 2 exp {ik [R(s, t) - R(s', if)]} iR^s', t') 
fc^iy^^exp^k [R{s,t) - R(s / ,t')]}^(s,t) 



kjkj\V(k)\ 2 exp {ik [R(s, t) - R(s', if)]} 
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I 6 (s,s';t,t') = (iRjis^iR^t) 



5R n (s,t)5R n (s',t') 



d 3 k 



kjki\V(k)\ 2 exp {ik[R(s,t) -~R(s',t'))} J. (B2) 
represent the selfconsistent averaging with the 



(2tt) 3 

The pointed brackets in eq 
Gaussian Hartree action. Taking this into account and using the generalized Wick 



theorem p2| , after straightforward algebra, we have 

1 r d 3 k ( k 2 "l 

h(s,s';t,t') = -G(s,s';t,t') J ±^k±\V{k)\* exp s'; t, t') j 

1 r d 3 k ( k 2 1 

I 2 (s,s';t,t f ) = -G(s,s f ;t',t) J — fc 4 |y(A;)| 2 exp \-jQ(a, a'; t,t')\ 

I 3 (s,s'-t,t') = ~G(s,s'-t,f) J ^k*\V(k)\ 2 exp^Q(s,s'-t,t')} 
= -h(s } s';t } t') 

h(s,s';t,t') = -~G(s,s>;t',t) J -0^k 4 \V(k)\ 2 exp {-yQ0>, s'; t, f)| 



t > t' 
t' > t 
t > t' 



t' > t 



I 5 (s,s';t,t') 



I s (s, s'-t,t') = 



-h(s,sf;t,t 
d 3 k 



k 2 



i) ^\V(k)\ 2 exp^- Y Q(s,s';t,t') 



(B3) 



where the last equation comes from the fact that the response function G(t, f) oc 
Q(t - t') and G(t', t) oc Q(t' - t). 



The 3-rd term in the exponent eq. (12) can be handled in the same way. The 
response function for the isotropic matrix has the form 



Pi(M) = -^P(M) 



(B4) 



where P(k, t) is the longitudinal part of the matrix response function. Then the 
Hartree approximation of the 3-rd term in the exponent ( |T2"D takes the form 

/ dsds / dt dt' {iR j (s,t)R j (s',t')J 1 (s,s';t,t')+iR j (s,t)R j (s,t)J 2 (s,s';t,t')} 

J0 J —oc J —oo ^ ' 
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where 



J 1 (s,s';t,t') = -J ^P|y(A;)| 2 P(k;t,t')exp|-^g( S)S ';t,0| 

r d 3 k ( k 2 1 

Ms,s';t,t') = J ^A; 2 |^)| 2 P(k;t,t')exp|- y g( S , S ';t,t')| 

= -J 1 (s,s';t,t') (B5) 



Taking into account eq. ( [BID with eq. ( [B2D and eq. ( |FJ5| ) leads to the Hartree-type 
approximation (|H>1). 
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FIGURES 



Renormalized static correlators (T=0.3) 
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FIG. 1. The static Rouse mode correlators for different chain lengths. The lines rep- 



resents the simple Rouse case. The temperatures are measured in units of the interaction 



potential with 7 = 0" = 1. 



29 



Critical temperature Tc(p/N) 




0.008 0.01 0.012 0.014 0.016 0.018 0.02 0.022 0.024 

p/N 

FIG. 2. The freezing temperatures of the Rouse mode correlators C(p) for different 
wavevectors p/N, where the temperatures are measured in units of the interaction potential 
with 7 = cr = l. The dashed line represents the freezing temperatures, when only the 



entropic contributions to the elastic susceptibility are taken into account. 
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Critical temperature for p=0,1 
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FIG. 3. The freezing temperatures for the p = and the p = 1 Rouse mode correlators 



C(p) of the test chain. 
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